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Abstract 


Experimental  and  theoretical  studies  of  the  optical  properties  of  large  core  step  index  (SI) 
plastic  optical  fibers  (POP)  and  graded  index  (GI)  POPs  are  reported.  A  set  of  criteria 
and  analyses  of  physical  parameters  are  developed  in  context  to  the  major  issues  of  POP 
applications  in  short-distance  communication  systems.  Analyses  are  presented  to  show 
how  the  measured  POP  optical  attenuation  affects  the  overall  performance  in  wavelength 
division  multiplexing  (WDM)  and  how  use  of  perfluorinated  polymers  can  overcome 
limitations  inherent  to  current  POP  materials.  Results  of  POP  optical  bandwidth 
measurements  by  direct  picosecond  time  domain  methods  are  reported  and  their 
relationship  to  refractive  index  profiles  theoretically  analyzed  by  the  WKB  and  finite- 
element  methods.  Two  high  resolution  optical  techniques  of  refracted  near-field  and 
transverse  interferometric  methods  are  presented  and  used  to  measure  the  index  profiles 
of  large  core  POPs.  Results  reveal  that  the  index  profile  of  current  GI  POP  is  not 
parabolic  ,  but  consists  of  two  markedly  different  regions.  Analysis  of  the  index  profile 
reveals  that  strong  mode  coupling  increases  the  GI  POP  bandwidth  from  its  profile- 
determined  value  of  0.43  GHz  to  its  measured  value  of  3.0  GHz  for  100m.  We  further 
study  mode  coupling  effect  inside  POPs  through  various  experimental  techniques, 
including  pulse  broadening  -  length  dependence  measurements  and  far-field  radiation 
pattern  measurements. 


Section  I.  Introduction 

As  the  role  of  organic  systems  in  photonics  increases,  considerable  attention  is 
centering  on  polymer  opto-electronic  devices  and  plastic  optical  fibers  (POFs),  especially 
in  fiber  link  applications  proposed  for  short  distance  (<500m)  communication  systems'. 
These  include  high  bandwidth  local  area  networks  (LANs),  data  links,  and  multi-node  bus 
networks^.  For  fiber  links,  the  key  point  is  that  installation  requires  many  junctions  and 
numerous  fiber  splices  and  connections,  and,  thus  ease  of  fiber  installation  and  handling 
become  important  considerations.  Typical  single  mode  glass  optical  fiber  (GOF)  has  a 
core  radius  of  5-10  p.m,  making  alignment  laborious  and  thus  installation  relatively  higher 
cost.  The  proposal  has  been  that  multi  mode  (MM)  POFs  can  satisfy  requirements  for 
quick-connect  fiber  links  for  short  distance  communications  systems.  The  key  advantage 
is  the  large  core  size  (~1  mm)  that  is  technically  and  economically  feasible  with  polymer 
materials.  Such  a  large  core  diameter  enables  quick  and  easy  connection  without  the  need 
for  precise  alignment  and  coupling  techniques.  Additionally,  large  core  POFs  are  easy  to 
handle  because  of  their  good  ductility  and  light  weight. 

Actual  selection  of  POF  fiber  link  technology  for  a  system,  however,  is 
determined  by  requirements  set  for  transmission  length  and  data  rate^.  To  date,  MM 
POFs  have  been  limited  to  relatively  low  data  rates  (<100  MHz- 100m)  and  short  lengths 
(<50m),  primarily  due  to  intermodal  dispersion  associated  with  the  step  index  (SI)  profile 
of  the  fiber  core.  In  GOFs,  the  well  established  technique  to  increase  bandwidth 
performance  is  to  grade  the  core  index  to  form  a  profile  approximate  to  a  parabola  such  as 
in  the  current  standard  62/125  MM  GOF  data  communication  fiber  links  that  can  operate 
at  greater  than  several  GHz-lOOm'*.  The  graded  index  (GI)  compensates  for  the 
intermodal  dispersion  and,  consequently,  increases  the  bandwidth  by  several  orders  of 
magnitude.  Co-polymerization  methods  have  been  proposed  for  realizing  large  core  MM 
GI  POFs^.  The  bandwidth  of  such  GI  POFs  has  been  independently  demonstrated  in  our 
laboratories  to  be  as  high  as  3  GHz- 100m  which  meets  a  main  requirement^. 

In  this  report,  we  summarize  the  results  of  experimental  and  theoretical  studies  of 
key  optical  properties  and  performance  levels  of  both  large  core  SI  POFs  and  GI  POFs 
based  on  polymethylmethacrylate  (PMMA).  A  primary  goal  is  to  provide  a  set  of  criteria 


and  analyses  of  major  issues  that  are  critical  to  utilizing  POFs  in  short-distance 
communication  systems.  The  first  sections  II-V  describe  the  early  phase  of  our  work  and 
introduce  major  issues  that  led  to  discoveries  concerning  mode  propagation  in  SI  POFs 
and  GI  POFs.  These  latter  observations  are  then  treated  separately  in  turn  in  Sections  VI 
and  Vn.  Sections  VI  and  Vn  contain  separate  references  in  addition.  In  section  n,  we 
first  analyze  how  measured  POF  transmission  properties  affect  important  issues  in 
wavelength  division  niultiplexing  (WDM)  and  discuss  how  perfluorinated  polymer 
materials  can  provide  a  path  to  overcome  inherent  limitations  of  current  POF  materials. 
Section  HI  contains  time-domain  bandwidth  measurements  of  POFs  and  theoretical 
analyses  of  the  POF  bandwidths  and  their  relationship  to  refractive  index  profiles. 
Section  IV  describes  two  different  high  resolution  optical  techniques  that  we  refined  to 
measure  large  core  POF  refractive  index  profiles.  Section  V  provides  a  summary  for 
Sections  II-IV. 

Section  11.  Wavelength  Division  Multiplexing  and  Optical  Attenuation 

There  exists  an  ever-growing  need  for  bandwidth  as  users  demand  more  and  more 
large  blocks  of  information.  The  bandwidth  required  per  user  has  recently  grown 
dramatically,  fueled  mainly  by  the  exploding  use  of  the  Internet  and  the  World  Wide 
Web.  To  handle  the  increasing  data  rates,  optical  fibers  are  becoming  the  links  of  choice 
in  current  communication  systems.  It  is  envisioned  that  optical  fiber  will  eventually 
displace  copper  wire  in  nearly  all  local  area  networks  (LANs)  and  will  even  reach  the 
home.  Therefore,  high  bandwidth  technologies  predicated  on  optical  fiber  have  become 
very  important  in  the  design  of  current  and  future  network  systems.  The  most  promising 
candidate  for  these  communication  systems  is  wavelength  division  multiplexing 
(WDM)l 

WDM  is  used  in  optical  communication  systems  to  effectively  increase  the 
transmission  bandwidth.  WDM  involves  multiplexing  a  number  of  channels  at  different 
wavelengths  onto  a  fiber  with  each  channel  operating  at  a  data  rate  within  the  dispersion- 
limited  bandwidth  of  the  fiber.  Consequently,  the  bandwidth  of  a  fiber  system  is 
effectively  increased  n-fold  for  an  n-channel  WDM  scheme.  Figure  4  schematically 


illustrates  the  concept  and  outlines  the  main  components  of  a  5-channel  WDM  system. 
The  diagram  shows  five  input  signals  at  different  wavelengths  written  from  five  separate 
laser  sources.  A  multiplexer  then  combines  the  signals  onto  a  single  optical  fiber.  For 
long  haul  systems,  the  signals  must  be  regenerated  or  amplified  before  reaching  their 
destination.  At  the  far  end,  the  five  signals  are  demultiplexed  and  directed  to  separate 
detectors  for  reading. 

WDM  is  important  for  both  point-to-point  communication  links  and  optical 
networks.  In  point-to-point  systems,  WDM  can  enhance  the  carrying  capacity  of  a  given 
fiber,  and  hence,  is  extremely  important  for  upgradability.  Additionally,  it  reduces  the 
number  of  fibers  required  to  support  a  particular  data  rate,  and  this  is  important  for 
situations  with  weight  and  volume  considerations.  Optical  networking  can  employ  a 
WDM  scheme  which  assigns  each  user  a  particular  wavelength.  For  example,  in  the  8- 
terminal  star  network  of  Figure  5,  each  user  transmits  one  wavelength  and  can  read  all  8 
wavelengths;  this  allows  a  user  to  communicate  with  any  other  user. 

The  factor  by  which  WDM  can  increase  the  bandwidth  of  a  fiber  system  is  only 
limited  by  the  number  of  channels  that  can  be  effectively  multiplexed.  Limitations  on  the 
number  of  channels  are  governed  by  the  wavelength  spacing  between  adjacent  channels 
and  the  total  wavelength  range  over  which  the  channels  can  be  operated.  The  minimum 
wavelength  spacing  is  determined  by  the  demultiplexing  detection  scheme,  the  source 
laser  linewidths,  and  nonlinear  optical  effects.  For  practical  WDM  systems  using 
relatively  inexpensive  components,  wavelength  spacings  of  one  to  several  nanometers  are 
achievable.  For  example,  demultiplexers  using  fixed  interference  filters  and  Fabry-Perot 
filters  to  isolate  wavelength  channels  have  resolutions  of  around  1  nm  and  distributed 
feed  back  laser  sources  can  have  linewidths  of  less  than  1  nm.  In  such  practical  systems, 
the  number  of  channels  is  determined  by  the  wavelength  range;  the  range  being  primarily 
limited  by  the  spectral  transmission  of  the  host  and,  in  long-haul  communications 
systems,  the  amplifier  gain  spectrum.  Consequently,  the  optical  attenuation  spectra  of 
candidate  fiber  media  are  of  critical  importance  to  WDM  system  performance. 


Optical  attenuation 

The  attenuation  spectrum  of  a  POF  can  be  experimentally  obtained  using  the 
standard  cut-back  technique.  The  attenuation  coefficient  Ot  in  dB  for  a  fiber  length  of  L 
is  given  by 

a,  =-Ylog,o(^).  (1) 

where  Pi  and  P2  are  the  input  and  output  powers,  respectively.  The  block  diagram  of 
Figure  1  shows  the  schematic  layout  of  the  experimental  setup®.  The  output  from  a  high 
power  xenon  lamp,  which  provides  an  incoherent  continuous  light  source,  passes  through 
a  computer-controlled  double-pass  monochromator.  After  being  modulated  by  a  chopper, 
the  light  from  the  monochromator  is  focused  onto  the  front  surface  of  the  POF  by  a 
convex -convex  lens  and  a  microscope  objective  lens.  The  transmission  from  the  fiber  is 
then  focused  onto  a  wide-area  silicon  photodiode  (PD).  Suitable  neutral  density  (ND) 
filters  are  placed  in  front  of  the  photodiode  to  avoid  saturation.  The  electrical  signal  from 
the  PD  is  then  sent  to  a  lock-in  amplifier  which  is  referenced  by  the  chopper.  The  output 
from  the  lock-in  amplifier  is  finally  fed  to  the  computer  for  automatic  data  acquisition. 

The  fiber  end  preparation  is  crucial  to  maintaining  the  accuracy  of  the 
measurement.  The  PDFs  in  our  measurements  were  polished  by  successively  using  four 
dry  lapping  films  of  decreasing  size  of  5,  3,  1,  and  0.3  |i.m,  respectively.  This  polishing 
technique  eliminates  end  surface  smearing  and  fiber  distortion  that  commonly  occur  with 
other  procedures  such  as  the  hot  plate  method. 

The  attenuation  spectra  of  a  typical  SI  POF  and  GI  POF  are  given  in  Figure  2. 
Since  both  the  SI  and  GI  POFs  are  PMMA-based,  there  are,  as  expected,  three  loss 
windows  located  at  570,  650,  and  780  nm  between  absorption  maxima  arising  from 
carbon-hydrogen  (CH)  vibrational  overtone  modes.  The  SI  POF  shows  excellent 
transmission,  having  a  loss  of  only  1 10  dB/km  at  the  principle  650  nm  window,  which  is 
very  close  to  the  theoretical  limit  of  106  dB/km.  The  GI  POF  also  possesses  satisfactory 
transmission  in  exhibiting  a  loss  of  158  dB/km  at  650  nm.  The  slight  increase  in  the  short 
wavelength  region  of  the  GI  POF  loss  spectrum  is  related  to  the  benzyl  benzoate  dopant 
used  in  forming  the  graded  index  profile. 


The  various  factors  contributing  to  optical  loss  in  POFs  can  be  divided  into 
intrinsic  and  extrinsic  loss’.  Intrinsic  loss  includes  vibrational  overtones  of  the  polymer 
materials,  electronic  transition  absorption,  and  Rayleigh  scattering.  Extrinsic  loss  includes 
absorption  due  to  impurities,  scattering  from  dust  and  microvoids,  and  imperfections  in 
fiber  parameters.  Intrinsic  loss  is  material  related  and  cannot  be  reduced  without  drastic 
changes  in  material  composition.  Extrinsic  loss  is  related  to  materials  processing  and 
fiber  fabrication,  and  thus  can  be  reduced  by  perfecting  each  procedure. 

Molecular  vibrations  of  aliphatic  hydrocarbons  are  the  dominant  intrinsic  loss 
factor  in  PMMA-based  POFs.  At  the  650  nm  window,  CH  absorption  contributes  about 
90  dB/km  to  the  total  loss.  Adding  in  electronic  absorption  caused  by  the  transition  of  the 
n  -K*  orbital  of  the  double  bond  within  the  ester  group  and  Rayleigh  scattering  caused  by 
density  fluctuations  of  the  polymer  matrix  brings  the  total  theoretical  loss  of  PMMA- 
based  POFs  to  106  dB/km  at  650  nm^. 

In  Ref.  3,  it  was  pointed  out  that  fluorinated  polymer  materials  provide  a 
promising  pathway  for  significantly  improving  transmission  performance  of  POFs  and 
here  we  show  how  this  is  important  to  meeting  critical  requirements  for  WDM.  Figure  3 
compares  the  theoretical  loss  for  a  perfluorinated  POF  with  the  above-measured  SI  and 
GI  PMMA  based  fibers.  The  primary  contributing  factors  to  the  perfluorinated  POF  loss 
are  wavelength-independent  structural  waveguide  imperfections  (10  dB/km)  and  Rayleigh 
scattering,  =95(568/  A)"*,  as  based  on  measurements  of  a  partially  fluorinated  core 
POF.  In  fact,  full  fluorination  would  reduce  the  polymer  index  of  refraction  and  turbidity 
and,  thus,  would  lead  to  an  even  lower  scattering  loss.  Losses  attributed  to  electronic 
absorption  are  negligible  as  are  absorptions  from  molecular  vibrations.  The  C-F 
vibrational  overtones  are  basically  nonexistant  in  the  visible  region  and  have  strengths  of 
much  less  than  1  dB/km  even  up  to  about  1500  nm®.  Consequently  the  total  loss  of  a 
perfluorinated  POF  should  be  less  than  25  dB/km  over  most  of  the  visible  spectrum  and 
approaches  10  dB/km  well  into  the  near  infrared. 


WDM  windows  of  operation:  Glass  and  plastic 

As  previously  discussed,  the  number  of  channels  in  a  WDM  system  is  limited  by 
the  widths  of  the  spectral  transmission  windows  of  the  fiber  medium.  In  the  case  of  silica 
glass  optical  fiber,  the  main  windows  of  operation  are  centered  at  the  1300  nm  and  the 
1550  nm  absorption  minima  (Figure  6).  There  is  little  chromatic  dispersion  around  1300 
nm  and,  therefore,  the  bandwidth  of  a  single  channel  operating  in  this  window  can  be 
very  high.  However,  the  advantages  of  the  1550  nm  window  are  that  it  can  support  more 
channels  with  its  broader  width  and  that  erbium  doped  fiber  amplifiers  can  be  used  to 
simultaneously  amplify  many  channels  in  long-haul  WDM  systems.  The  wavelength 
range  for  WDM  is  limited  by  the  width  of  the  transmission  window  because  significant 
differences  in  attenuation  between  channels  will  result  in  large  differences  in  signal 
intensities  at  the  detectors.  For  example,  if  one  channel  experiences  an  attenuation  of  0.3 
dB/km  while  another  channel  suffers  a  loss  of  0.5  dB/km,  then  after  100  km  the  two 
signal  intensities  would  differ  by  20  dB  which  is  clearly  not  acceptable. 

In  fact,  intensity  differences  between  channels  of  less  than  3  dB  are  necessary  for 
similar  bit  error  rates  to  be  achievable,  therefore,  we  define  the  width  of  a  particular 
transmission  window  to  be  the  range  of  wavelengths  for  which  there  would  be  less  than  a 
3  dB  intensity  difference  between  signals  at  the  window  boundaries,  ^boundary,  and  a  signal 
at  the  window  absorption  minimum,  Xmi„,  over  a  distance,  L,  for  which  the  latter  suffers 
an  attenuation  of  20  dB.  In  other  words 


A(.A^)Y  L=3dB 


(2) 


where 


L  =  20dBI 


(3) 


and  A(k)  is  the  attenuation  of  the  fiber  in  dB  at  the  wavelength,  X.  Accordingly,  the  1550 
nm  window,  which  has  a  width  of  180  nm,  is  preferable  to  the  1 10  nm  wide  1300  nm 
window  as  it  would  allow  nearly  twice  as  many  WDM  channels  (Table  1). 

In  contrast  to  silica  glass  fibers,  PMMA-based  PDFs  have  a  very  narrow 
transmission  window  (Figure  2)  of  only  approximately  10  nm  at  the  proposed  operating 
wavelength  of  650  nm  (Table  1).  Therefore,  the  number  of  possible  channels  for  WDM 
is  significantly  decreased  in  POF  systems.  For  example,  in  a  WDM  system  with  channels 
spaced  4  nm  apart,  silica  glass  fibers  operating  within  the  1550  nm  window  could 
support  up  to  45  different  channels  while  650  nm-centered  PMMA  PDFs  could  handle 
only  3.  Such  a  dramatic  difference  in  performance  exposes  the  inherent  limitations  on 
WDM  applications  based  around  the  650  nm  operating  wavelength. 

We,  therefore,  conclude  that  if  WDM  is  to  be  viable  in  POF  systems,  either  the 
operating  wavelength  in  PMMA-based  POFs  should  be  moved  to  the  wider  transmission 
windows  below  600  nm,  or  POFs  should  be  based  on  other  materials  such  as  fluorinated 
polymers  that  have  lower  losses  and  broader  windows.  Analyzing  the  lower  wavelength 
windows  (Figure  2)  for  both  step-index,  SI,  and  graded-index,  GI,  PMMA-based  POFs, 
we  can  see  that  there  is  not  a  significant  improvement  in  the  number  of  WDM  channels 
(Table  2).  The  SI  fibers  can  accommodate  up  to  a  maximum  of  10  channels  in  the 
transmission  window  centered  at  520  nm.  However,  this  window  is  completely 
eliminated  in  GI  POFs  due  to  large  absorption  by  the  benzyl  benzoate  dopant  used  in 
forming  the  index  gradient.  Consequently,  the  570  nm  window  is  the  widest,  allowing  up 
to  8  distinct  channels. 

Fluorinated  fibers  would  provide  broad  a  transmission  window  extending  over  the 
entire  visible  region  (Figure  3)  for  enhanced  WDM  operation.  In  addition,  fluorinated 
POFs’  would  possess  significantly  lower  loss  in  the  both  the  visible  and  infrared  regions 
and,  consequently,  would  allow  operation  over  a  wide  spectral  range  including  the  850, 
1300,  and  1550  nm  communication  wavelengths.  These  fluorinated  polymer  features  are 
more  desirable  than  the  lower  wavelength  PMMA  windows  as  the  lack  of  high 
performance  transceivers  below  600  nm  would  preclude  such  operating  wavelengths. 


Section  III.  Bandwidth 
Time  Domain  Measurement 

The  bandwidth  performance  of  an  optical  fiber  can  be  determined  by  measuring 
either  the  impulse  response  in  the  time  domain,  or  the  spectral  response  in  the  frequency 
domain.  In  our  experiments,  we  measured  the  impulse  response  of  a  GI  POP  in  the  time 
domain  using  a  high  repetition  rate,  picosecond  laser  diode.  Picosecond  (ps)  light  sources 
are  ideal  for  the  measurement  since  nanosecond  pulses  do  not  have  the  required 
resolution.  Also,  the  use  of  femtosecond  pulses  might  inadvertently  introduce  additional 
broadening  due  to  self-phase  modulation.  In  our  current  setup,  with  X\  and  X2  the  full 
widths  at  half  maximum  (FHWM)  of  the  light  source  and  output  signals  of  the  fiber, 
respectively,  the  impulse  response  at  FWHM  of  the  fiber  is  given  by 

.  (4) 


Furthermore,  if  all  pulses  involved  are  assumed  Gaussian,  the  fiber  bandwidth  is 


B1T  = 


0.44 

- L, 


(5) 


where  L  is  the  fiber  length.  In  general,  however,  fast  Fourier  transforms  are  required  to 
obtain  accurate  bandwidth  measurements. 

Our  time  domain  bandwidth  measurement  setup  is  shown  in  Figure  7.  The  light 
source  is  a  temperature  stabilized  InGaAlP  laser  diode  capable  of  producing  45  ps  pulses 
at  660  nm  wavelength.  The  repetition  rate  of  the  laser  diode  can  be  as  high  as  2  MHz. 
The  laser  diode  output  is  first  collimated  by  a  GRIN  lens  and  spatially  filtered.  It  is  then 
collimated  by  a  15  cm  focal  length  lens  and  focused  onto  the  POF  by  a  microscope 
objective  lens  (lOX,  0.25  NA).  The  light  exiting  the  other  end  of  the  fiber  is  detected  by  a 
sampling  optical  oscilloscope.  The  time  delay  of  the  trigger  signal  for  the  optical 
oscilloscope  can  be  adjusted  by  the  delay  box  as  well  as  the  diode  controller.  The  main 
part  of  the  optical  oscilloscope  consists  of  a  sampling  streak  tube  which  converts  the  light 


incident  on  its  photocathode  into  photoelectrons.  The  photoelectrons  are  sampled  and 
directed  to  a  phosphor  screen,  where  they  are  converted  to  light.  This  light  is  then 
detected  by  a  photomultiplier  tube  (PMT).  The  data  stored  in  the  oscilloscope  is 
transferred  to  the  computer  for  further  data  analysis.  The  time  resolution  of  the  system  is 
better  than  10  ps. 

In  the  Gaussian  approximation,  the  effective  NA  of  the  light  beam  incident  onto 
the  fiber  can  be  determined  by 


NA  = 


mfi 


(6) 


where  X  is  the  wavelength,/]  (15  cm)  is  the  focal  length  of  the  collimating  lens,  a  (25 
pm)  is  the  pinhole  radius  for  the  spatial  filter,  and  /2  (14.8  mm)  is  the  focal  length  of  the 
microscope  objective  lens.  The  effective  NA  is  thus  calculated  to  be  0.09.  Since  all  of  the 
POP  samples  that  we  have  tested  have  NAs  of  near  0.2,  they  are  underfilled  by  the  input 
light.  An  underfilled  launch  condition  tends  to  produce  higher  bandwidth  than  an 
overfilled  condition  because  only  a  portion  of  the  modes  are  excited,  resulting  in  a 
smaller  total  dispersion. 

The  dependence  of  the  measured  bandwidth  on  the  launch  condition  is  solved  by 
introducing  a  mode  scrambler  immediately  after  the  launching  end  of  the  sample  fiber. 
The  function  of  the  mode  scrambler  is  to  induce  mode  couplings  among  all  groups  of 
modes  by  introducing  periodic  bending  in  the  fiber  so  that  an  equilibrium  mode 
distribution  is  reached.  An  illustration  of  the  mode  scrambler  that  we  designed  for  use 
with  POPs  is  shown  in  Pigure  8.  Our  mode  scrambler  consists  of  two  corrugated  surfaces 
between  which  the  fiber  is  placed.  The  optimal  spatial  period  is  determined  by  the 
expression 


where  d  is  the  fiber  core  diameter  and  A  is  the  fractional  refractive  index  difference 
between  the  core  center  and  cladding.  The  typical  values  of  d  and  A  for  GI  POP  are  500 
|xm  and  0.01,  respectively,  resulting  in  a  length  of  1 1  mm  for  the  spatial  period  A. 

To  verify  the  effectiveness  of  the  mode  scrambler,  we  performed  the  following 
experiment  with  and  without  the  mode  scrambler.  For  a  98  m  Mitsubishi  GI  POP  with  a 
NA  of  0.21  and  an  ttj  of  158  dB/km,  we  measured  the  pulse  width  of  the  output  while 
varying  the  incidence  angle  of  the  launch  beam.  The  results  are  shown  in  Figure  9. 
Without  the  mode  scrambler  in  place,  the  pulse  width  of  the  output  pulse  varies 
significantly  with  the  incident  angle  because  at  different  incident  angles,  different  .groups 
of  modes  are  excited  which  exhibit  different  dispersions.  With  the  mode  scrambler,  the 
variation  in  the  pulse  width  is  much  smaller,  indicating  that  an  equilibrium  mode 
distribution  is  achieved.  It  is,  therefore,  necessary  to  use  the  mode  scrambler  in  the 
bandwidth  measurement  to  achieve  consistent  results. 

We  have  measured  both  SI  POP  and  GI  POP  samples  with  our  time-domain  setup. 
Figure  10  shows  the  input  and  output  pulses  for  a  Mitsubishi  SI  POP  with  a  NA  of  0.47 
and  ttT  of  1 10  dB/km.  The  bandwidth  of  this  particular  fiber  is  found  to  be  80  MHz  for  a 
100  meter  length.  Figure  11  shows  the  input  and  output  pulses  for  a  98  m  Mitsubishi 
Rayon  GI  POP  with  a  NA  of  0.21  and  aj  of  158  dB/km.  The  bandwidth  is  determined  to 
be  3.0  GHz  for  100  m.  All  measurements  are  performed  with  the  mode  scrambler  applied 
to  the  fiber  sample. 

Bandwidth  Analysis  and  Refractive  Index  Profile 

There  exist  several  different  techniques  to  analyze  the  bandwidth  properties  and 
performance  of  a  MM  optical  fiber.  These  include  the  WKB  approximation,  Rayleigh- 
Ritz,  power-series  expansion,  and  finite-element  methods.  For  an  SI  fiber,  one  can  obtain 
results  analytically  using  the  following  equation'®  which  is  derived  by  requiring  that 
waves  are  continuous  at  the  core-cladding  boundary: 


Uj'y{u)-W  =0, 

K^{w) 


(8) 


where  and  are  Bessel  functions  of  the  first  kind  and  modified  Bessel  functions  of  the 
first  kind,  respectively,  v  is  the  azimuthal  number.  The  normalized  propagation  constant  u 
is  related  to  the  ko,  the  wave  vector  in  vacuum,  rif,  the  core  refractive  index,  (3,  the 
propagation  constant,  and  a,  the  radius  of  the  fiber  core  by  the  following 

u  =  ^ikln]-P^)a.  ■  (9) 


w  is  defined  by 


w 


=  ylV^-U^ 


(10) 


where  V  is  the  normalized  frequency  number. 


V  =  ^(k^nj  -klnl)a  =  . 


(11) 


A  is  the  relative  index  difference  between  the  core  refractive  index  rif  and  the  cladding 
refractive  index  tic,  and  is  given  by  A  =  {nl  -n])  I (2nj ) . 


The  group  delay  time  x  can  be  calculated  through 


10 


1-2A0 


T  =  - 


(r 


(12) 


where 


(13) 


0  = 


y/{W)  - 1 


+  w(yv) 


In  the  above  equation,  (l)Bf w)  and  V|t( w)  are  determined  by 


Ij/iw)  = 


Kliw) 


(14) 

(15) 


The  impulse  response  function  can  then  be  constructed  by  counting  the  number  of  modes 
in  each  time  interval,  taking  into  consideration  that  the  degeneracy  for  modes  with  v  =  0 
is  2  and  v  0  is  4.  Figure  12  shows  the  resultant  impulse  response  function  for  an  SI 
fiber  of  V  =  2000.0  and  A  =  0.01,  assuming  all  modes  are  equally  excited.  The  maximum 
group  delay  time  is  An/c,  which  is  the  group  delay  time  difference  between  the  lowest 
mode  traveling  with  velocity  cMf  and  the  highest  mode  traveling  with  velocity  cMc.  The 
main  feature  of  the  impulse  response  function  is  that  it  is  similar  to  a  square  wave,  with  a 
slight  downward  slope  as  the  group  delay  time  increases,  a  feature  we  will  explain 
shortly.  The  spike  near  the  cut-off  is  due  to  the  high-order  modes. 

WKB  approximation 

The  semi-classical  WKB  approximation  is  a  widely  used  technique  to  study 
optical  properties  of  fibers.  It  is  especially  appropriate  for  MM  POFs,  where  the  number 
of  modes  is  on  the  order  of  hundreds  of  thousands  so  that  the  semi-classical  result  is  quite 
accurate.  Under  the  weakly  guiding*’  and  slowly  varying  approximations,  the  following 
scalar  wave  equation  can  be  obtained. 


( 

\p- 

LA. 

r  dr 


^  =  0, 


(16) 


where  v  is  the  azimuthal  number  and  R(r)  is  the  radial  part  of  the  electric  field.  The 
impulse  response  function  h(x)  can  be  calculated  using  the  following  equations: 


^  ir)ko  -  J3^)rdr 

(17) 

1  dm  /  dkn 

^  c  dm  1  dp 

(18) 

dm 

h(,t)  = 

(19) 

where  m  is  the  number  of  modes  above  the  propagation  constant,  P,  h  is  the  wave 
number,  rj  is  determined  by  n{r2)  =  P  /  /^o,  and  X  is  the  group  delay  time,  and  p{}n)  is  the 
modal  power  distribution.  Later,  we  will  implement  this  technique  to  analyze  the  impulse 
response  for  a  GI  POP  from  its  measured  refractive  index  profile.  For  an  SI  POP,  it  can 
be  shown  that  the  WKB  impulse  response  function,  assuming  a  uniform  modal  power 
distribuition,  p{m)  =  1,  is 


/j(t)  oci/  for  I  c<T<n^f  I  (n^c) . 


(20) 


The  group  delay  time  can  be  written  as  T  =  t/+  8t  where  Xf=nf  I  c  is  the  delay  time  of  the 
lowest  guided  mode  and  6t  is  the  delay  time  relative  to  Xf.  Since  6x  «  X/,  h(x)  satisfies 


h{T)oc\-3ST/ Tf,  0<St<—A. 

’  c 


(21) 


The  above  equation  demonstrates  that  the  impulse  response  function  for  an  SI  fiber 
behaves  as  a  square  wave  with  a  slight  negative  slope  in  its  height,  in  excellent  agreement 
with  our  previous  analytical  results. 


Finite  element  method 

We  have  also  developed  a  finite  element  method  for  POF  as  adapted  from  that 
used  by  Okamoto*^  for  GOF,  in  which  the  region  between  r  =  0  and  r  =  a  is  divided  into 
N  elements,  and  the  values  of  the  radial  part  of  the  electric  field  R(r)  at  r  =  r,  are 
expressed  as 


R.=R(r^),  i  =  0,1,2,..., N .  (22) 

We  can  then  construct  i?(r)  using  a  combination  of  linear  functions  to  produce  a 
tridiagonal  matrix  (see  Appendix  A).  The  eigenvalues  of  the  matrix  correspond  to  u, 
which  is  directly  related  to  the  propagation  constant  P,  and  the  eigenvectors  represent  the 
radial  distribution  of  the  electric  field.  The  group  delay  time  is  straightforwardly 
calculated  from  these  quantities.  We  have  carried  out  finite  element  calculations  on  fibers 
with  various  refractive  index  profiles.  Figure  13  shows  the  ratio  of  the  group  delay  time 
for  the  stair  case  profile  (top)  to  that  of  the  step  index  profile,  as  a  function  of  the  number 
of  steps  (bottom).  Such  stair  case  profiles  have  important  practical  implications  because 
they  can  be  easily  fabricated.  The  refractive  index  nfr)  is  determined  by 


for;jfl  <;ir<  1, 

(23) 

nc 

for2r>  1, 

where 

X-r  la  (24) 

An  interesting  feature  of  Figure  13  is  that  the  ratio  begins  to  saturate  after  only  about 
seven  steps. 

Using  the  finite  element  method,  we  are  also  able  to  calculate  the  radial  part  of  the 
electric  field  R{f)  for  any  given  refractive  index  profile.  Figure  14  shows  Rij)  of  a 
parabolic  GI  fiber  for  three  different  guided  modes,  all  of  which  are  characterized  by  (v, 
/),  where  v  is  the  azimuthal  number  and  I  is  the  sequential  number  in  the  v  guided-mode 


group.  The  number  of  elements  used  in  the  calculation  is  N  =  1000,  the  normalized 
frequency  number  V  =  400.0,  and  A  =  0.01.  These  graphs  clearly  show  that  for  low-lying 
guided  modes,  the  light  is  confined  to  a  certain  region  within  the  fiber,  consistent  with 
findings  from  ray  theory.  High-lying  guided-modes,  however,  are  more  oscillatory  and, 
therefore,  cannot  be  adequately  described  by  ray  theory. 

Section  IV.  Measurements  of  Refractive  Index  Profiles  of  Plastic  Optical  Fibers 
Refracted  Near-Field  Measurement 

There  are  several  different  methods  to  measure  the  refractive  index  profile  of  an 
optical  fiber.  A  commonly  used  one  is  the  refracted  near  field  technique*^  ''*  in  which  one 
measures  the  intensity  of  refracted  light  escaping  from  the  fiber.  Figure  15  shows  the 
schematic  experimental  setup‘s.  The  650  nm  CW  laser  light  is  spatially  filtered  with  a 
lOX  microscope  objective  lens  and  a  pinhole.  The  beam  is  then  converted  into  circularly 
polarized  light  by  the  combination  of  a  polarizer  and  a  quarter-wave  plate.  This 
minimizes  the  angular  dependence  of  the  reflectivity  at  the  glass-air  interface.  The  CW 
light  is  then  focused  onto  the  fiber  front  end  by  a  25X  microscope  objective  lens  that  has 
a  numerical  aperture  (NA  =  0.50)  much  greater  than  that  of  the  fiber.  The  refracted  light 
is  collected  onto  a  silicon  photodiode  by  a  Fresnel  lens.  Leaky  rays’^’’'*  are  prevented 
from  reaching  the  detector  by  the  blocking  disk.  The  fiber  is  immersed  in  an  index 
matching  liquid  contained  in  a  sample  cell  mounted  on  a  3-D  piezoelectric  stage.  The 
front  window  of  the  sample  cell  is  a  200  |xm  thick  microscope  slide  coverglass.  The  back 
window  is  a  plexiglass  plate  having  a  tiny  hole  in  the  center  through  which  the  fiber  is 
inserted.  Light  from  a  He-Ne  laser  is  injected  into  the  back  end  of  the  fiber  and,  in 
conjunction  with  an  eyepiece,  is  used  to  align  and  focus  the  650  nm  beam  onto  the  fiber’s 
front  endface.  By  focusing  the  light  onto  the  fiber  endface  and  then  scanning  the  endface 
across  the  beam,  the  local  refractive  index  can  be  measured  through  the  following 
relation 

P  ^  (sin'  max  -  sin'  +nl  -n\r)) 


(25) 


where  P  is  the  power  detected,  Gamin  is  the  angle  defined  by  the  blocking  disk,  G/mn^  is 
directly  related  to  the  numerical  aperture  of  the  input  light,  nt  is  the  refractive  index  of 
the  matching  liquid,  and  n(r)  is  the  local  refractive  index  of  the  fiber.  When  tit  ~  n{r),  the 
detected  power  P  is  a  linear  function  of  n{r). 

The  blocking  disk  is  crucial  to  obtaining  an  accurate  measurement.  For  this 
purpose,  we  have  constructed  a  blocking  disk  which  is  made  of  aluminum  and  attached  to 
a  large  thin  piece  of  plexiglass,  as  shown  in  Figure  16.  A  fiber  holder  is  placed  through 
the  hollow  bolt  to  feed  the  fiber  to  the  sample  cell.  The  diameter  of  the  blocking  disk  is 
two  inches.  The  large  diameter  enables  us  to  place  the  blocking  disk  further  from  the 
sample  cell  to  minimize  image  aberration. 

We  can  determine  the  absolute  refractive  index  through  two  calibration  points:  the 
known  refractive  indices  of  the  matching  liquid  (1.501)  and  the  PMMA  cladding  (1.491). 
The  measured  refractive  index  profile  is  plotted  in  Figure  17.  The  distinguishing  feature 
of  the  measured  profile  is  that  it  is  not  parabolic  as  one  would  expect  for  a  true  G1  fiber. 
The  central  part  of  the  profile,  near  to  the  core  center,  is  closely  parabolic,  but  the  outer 
region  connects  to  the  cladding  more  steeply  than  would  a  parabola. 

We  can  fit  the  measured  index  profile  to  a  polynomial  profile  given  by 


«(r)  = 


n. 


-  dn^^ir  !  a)^  -  !  aY  forr<a^, 

-  Sn2oir  I  a)^  -  Sn2i(r  /  a)^  fora^<r<a, 

for  r>  a. 


(26) 


We  have  found  that  two  polynomials,  corresponding  to  the  central-core  {r<ai)  and  core¬ 
cladding  {uj  <r<a)  regions,  are  necessary  to  properly  fit  the  measured  profile.  Requiring 
the  two  polynomials  to  connect  smoothly  at  r  =  ay,  the  fitted  parameters  are  =  1 .5 1 14, 

5n,o  =  0.0102,  5/Jii  =  0.0038,  =  1.5278,  5n20  =  0.0173,  bnzi  =  0.0197,  ric  =  1.4910,  a, 

=  235.0  |im,  and  a  =  260  |Xm.  The  fact  that  5^21  is  much  greater  than  dnn  demonstrates 
that  in  the  core-cladding  region,  the  refractive  index  drops  precipitously  to  the  cladding. 
The  polynomial  fit  is  also  shown  in  Figure  17.  Initially,  we  attempted  to  fit  the  measured 


data  to  other  profiles,  including  an  alpha  profile  and  even  a  single  sixth-order  polynomial, 
and  found  that  they  are  not  able  to  map  the  features  of  the  measured  data  throughout  the 
entire  core  region.  Figure  18  shows  the  measured  profile  along  with  an  alpha  fit  with  the 
exponent  of  a  =  3.67.  It  is  obvious  that  the  alpha  profile  does  not  agree  with  the  data  in 
the  central  core  region. 

In  addition  to  intermodal  dispersion,  material  dispersion  can  contribute  to  pulse 
broadening  in  a  MM  optical  fiber‘^.  The  bandwidth  of  a  fiber  with  a  near-parabolic 
refractive  index  profile  (a  =  2.0)  can  be  greatly  affected  by  material  dispersion. 
However,  we  find  that  the  for  alpha  profiles  with  a  >  3.0  and  source  spectral  widths  less 
than  5  nm,  the  effect  of  material  dispersion  is  negligible.  This  fact  is  independently 
evidenced  in  Fig.  7  of  Ref.  17. 

The  calculated  impulse  response  function  using  the  WKB  approximation  for  a  100 
m  GI  POF  having  the  measured  profile  is  shown  in  Figure  19.  For  our  calculation,  the 
modal  power  distribution  of  equation  19  is  taken  to  be'^ 


p(m)  = 


fl -(m/MJ 

lo 


for  m  < 
for  m  > 


(27) 


where  Me  represents  the  maximum  number  of  modes  that  can  propagate  in  the  fiber. 
Here,  Me  is  set  to  0.95M,  where  M  is  the  total  number  of  guided  modes.  The  impulse 
response  of  Figure  19  has  two  dissimilar  parts.  The  first  part  is  a  sharp  and  intense,  and  is 
the  result  of  the  propagation  of  modes  confined  to  the  central  core  region.  The  second 
part  is  very  broad  with  much  lower  intensity,  representing  the  contributions  from  the 
remaining  modes  which  travel  with  large  dispersion  in  both  the  central-core  and  core¬ 
cladding  regions.  The  Fourier  transform  of  the  impulse  response  function  yields  the 
bandwidth,  which  is  the  frequency  at  which  the  transfer  function  (see  Figure  20)  drops  to 
half  of  its  peak  value.  The  bandwidth  of  a  100  m  length  of  current  GI  POF  is  determined 
to  be  0.45  GHz  in  this  manner. 

For  comparision,  direct  measurements  of  fiber  bandwidths  using  the  time-domain 
method  yield  a  bandwidth  of  3.0  GHz  for  100  m,  which  is  much  higher  than  the 


calculated  bandwidth.  We  ascribe  the  difference  between  the  calculated  intrinsic  and 
measured  bandwidths  to  mode  coupling  effects.  After  equilibrium  mode  distribution 
(EMD)  is  established  through  mode  coupling,  the  pulse  broadening  decreases  from  being 
linearly  proportional  to  the  fiber  length  to  varying  as  the  square  root  of  length,  effectively 
increasing  the  bandwidth'^.  We  have  observed  mode  coupling  effects  in  POFs'^,  with  the 
result,  for  example,  that  the  EMD  condition  is  achieved  in  SI  PDFs  between  15  and  20  m. 
Preliminary  measurements  indicate  that  the  coupling  length  in  GI  POFs  is  even  shorter. 
Furthermore,  it  can  be  shown  that 


where  BWx  and  5W2  are  the  bandwidths  with  and  without  the  mode  coupling  effect, 
respectively,  L  is  the  fiber  length,  and  Lc  is  the  coupling  length.  Applying  equation  28  to 
our  experimental  and  calculated  bandwidths,  the  coupling  length  is  determined  to  be  2.2 
m,  consistent  with  our  preliminary  results. 

Transverse  interferometric  measurement 

A  second  measurement  was  developed  to  corroborate  the  GI  POF  index  of  refraction 
profile  obtained  using  the  near-field  refracted  ray  technique.  To  this  end,  a  transverse 
interferometric  method  (TIM)  was  adapted  to  measure  POF  refractive  index  profiles. 
Previously^”,  TIM  had  been  applied  to  single-mode  and  MM  glass  optical  fibers  having 
maximum  outer  diameters  of  125  microns.  Several  modifications  were  necessary  to 
allow  the  measurement  to  handle  thicker  POFs  up  to  and  exceeding  1  mm  in  diameter. 

The  transverse  interferometric  method  is  implemented  using  a  Leitz  transmitted- 
light  interference  microscope.  As  depicted  in  Figure  21,  a  beam  of  light  from  a  tungsten 
lamp  is  split  into  two  arms.  One  arm  serves  as  the  reference  while  the  other  passes 
through  the  fiber  sample,  transverse  to  the  fiber  axis.  The  two  beams  are  then 
recombined  and  interfere  with  each  other.  In  the  absence  of  a  fiber  sample,  a  uniform 
bright  output  will  result  if  the  path  lengths  of  the  two  arms  are  identical.  However, 


slightly  tilting  one  arm  with  respect  to  the  other  causes  a  pattern  of  evenly-spaced, 
parallel  fringes  to  appear  at  the  output.  The  displacement  of  these  fringes,  upon  insertion 
of  a  fiber  sample,  can  then  be  used  to  derive  the  radial  dependence  of  the  fiber  index  of 
refraction  profile.  The  output  is  viewed  through  the  microscope  eyepieces,  and  is 
simultaneously  recorded  by  a  CCD  camera  mounted  on  top  of  the  interference 
microscope.  The  CCD  camera  images  are  displayed  on  a  video  monitor.  In  order  to 
analyze  the  fringe  displacements,  the  camera  is  interfaced  to  a  computer  via  a  video  board 
that  captures  the  camera  output  and  converts  it  to  digital  images.  The  stored  fringe 
displacements  are  then  traced  and  used  to  calculate  the  fiber  index  of  refraction  profile. 

Figure  22  schematically  illustrates  the  interference  fringe  patterns  before  and  after 
insertion  of  the  fiber  sample.  The  evenly-spaced,  parallel  fringes  are  displaced  by  the 
presence  of  the  fiber.  The  index  of  refraction  differences  experienced  by  rays  passing 
through  the  fiber  cause  changes  in  the  effective  path  lengths  traveled  by  those  rays. 
Consequently,  a  given  ray  suffers  an  additional  phase  shift. 


A(1)  =  ^|J  [«(/)- 


(29) 


where  k=2KlX  is  the  wavenumber  corresponding  to  the  wavelength  X,  n{l)  is  the  refractive 
index  variation  along  the  path  (Figure  23),  /,  of  the  ray,  and  hl  is  the  refractive  index  of 
the  index  matching  liquid.  The  measured  fringe  shift,  S,  corresponding  to  a  phase  shift  of 
At])  follows  from 


D  2% 


(30) 


where  D  is  the  distance  between  fringes.  Extracting  S  and  D  from  the  electronically 
captured  fringe  pattern  allows  numerical  determination^ of  the  desired  refractive  index 
profile,  n(l)-ni. 

The  two  arms  of  the  microscope  stage  are  displayed  in  Figure  24.  On  the  sample 
side,  the  fiber  is  placed  on  a  quartz  plate  and  immersed  in  index  matching  fluid  that 


matches  the  fiber  cladding  index.  Sufficient  matching  liquid  is  used  so  as  to  completely 
cover  the  lens  of  the  oil-immersion  objective.  On  the  reference  side  the  objective  lens  is 
immersed  in  the  same  index  matching  fluid  over  an  identical  quartz  plate.  Whereas  lOOX 
microscope  objectives  are  used  for  glass  fibers,  20X  lenses  are  implemented  to 
accommodate  a  larger  field  of  view  as  necessary  for  the  thicker  POFs.  The  other  major 
difference  between  the  plastic  and  glass  measurements  is  that  the  fringe  displacements  are 
significantly  larger  for  POFs  owing  to  the  much  longer  path  lengths  traveled  through  the 
fiber  by  the  sample  beams. 

Experimental  results  for  the  GI  POF  index  of  refraction  profile  are  displayed  in 
Figure  25  and  compared  to  the  near  field  refracted  ray  measurement.  The  two  profiles  are 
nearly  identical,  confirming  the  accuracy  of  both  techniques.  Additionally,  the  theoretical 
bandwidth  calculated  for  this  TIM  profile  is  0.38  GHz  -  100m,  corroborating  the  0.45 
GHz  -  100  m  bandwidth  results  found  via  the  near  field  refracted  ray  measurement. 

Section  V.  Conclusion 

We  have  conducted  experimental  and  theoretical  studies  of  key  optical  properties 
and  performance  levels  of  both  large  core,  PMMA-based  SI  and  GI  POFs.  Additionally 
we  have  provided  a  set  of  criteria  and  analyses  of  major  issues  that  are  critical  to  utilizing 
POFs  in  short-distance  communication  systems.  In  particular,  we  have  described  how 
optical  attenuation  limits  fiber  bandwidth  and  transmission  length.  Our  measurements 
reveal  that  the  loss  of  current  PMMA-based  SI  POFs  is  approaching  the  theoretical  limit. 
Additionally,  absorption  due  to  electronic  excitations  of  dopant  molecules  utilized  in  GI 
POFs  can  increase  loss  and  even  eliminate  transmission  windows.  We  have  presented  an 
analysis  relating  the  optical  attenuation  spectrum  to  the  overall  performance  in  WDM 
systems.  We  conclude  that  hydrocarbon-based  polymers  will  not  qualify  for  such 
applications,  but  rather  fluorinated  materials  may  be  necessary.  We  have  performed  time- 
domain  bandwidth  measurements  of  both  SI  and  GI  POFs  and  have  related  the  results  to 
theoretical  analyses  of  their  refractive  index  profiles  using  WKB  and  finite-element 
methods.  Further,  we  have  developed  the  refracted  near-field  and  transverse 
interferometric  methods  for  measuring  the  refractive  index  profiles  of  large  core  POFs. 


Results  reveal  that  the  index  profile  of  current  GI  POF  is  not  parabolic  in  shape  which 
significantly  limits  the  bandwidth  compared  to  that  of  a  parabolic  GI  profile. 


Appendix  A 

In  the  bandwidth  analysis  using  the  finite-element  technique'^,  a  tridiagonal, 
symmetric  matrix  5,y  (A+1  by  A^+1,  where  N  is  the  number  of  elements)  is  constructed 
whose  eigenvalues  correspond  to  the  normalized  propagation  constant  u  and  whose 
eigenvectors  correspond  to  the  radial  part  of  the  electric  field  /?(r).  The  non-zero 
elements  of  the  symmetric  matrix  Sy  are 
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and  u,  V,  and  w  are  defined  in  Equations  (9),  (11),  and  (10),  respectively,  rif  is  the 
refractive  index  at  the  center  of  the  fiber  core.  A  is  the  relative  refractive  index  difference 
between  the  core  center  and  cladding,  and  v  is  the  azimuthal  number. 


J  2  • 

The  results  we  obtained  above  are  slightly  different  from  Okamoto’s  ,  especially 
for  5oo  and  5oi.  In  addition,  we  also  derived  the  correct  expression  for  Su,  in  which  the 
coefficient  associated  with  qo  should  be  3  instead  of  30. 

The  group  delay  time  x  can  in  general  be  calculated  by  the  following 

where  the  function  f(ko,  P)  =  0  determines  the  (ko,  p)  dispersion.  In  some  cases,  the 
dispersion  function  is  in  the  form  of  f(u,  v)  =  0,  in  which  case  the  group  delay  time  x  can 
be  determined  using  the  following; 

^  ^  du  ^  dv 

^  ^  du  ^  dv  ^  Bu 

since  v  is  independent  of  p.  We  therefore  derive  the  following 
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The  above  equation  is  particularly  useful  in  the  power-series  expansion  method^'^  for  a- 
profiles.  We  found  that  the  corresponding  equation  in  Ref.  23  was  incorrect.  The  relative 
group  delay  time  T  is  then 
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T  is  usually  small  since  A  «  1 . 
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Table  1.  Transmission  windows  for  silica  glass  and  PMMA  plastic  optical  fiber. 


Optical  Fiber 

Center 

Operating 

Wavelength 

Width  of 
Window 

Number  of  WDM 
Channels  (4  nm  spacing) 

Silica  Glass 

1550  nm 

180  nm 

45 

Silica  Glass 

1300  nm 

llOnm 

28 

PMMA  Plastic 

650  nm 

10  nm 

3 

Table  2.  Transmission  windows  for  SI  and  GI PMMA  plastic  optical  fiber. 


Optical  Fiber 

Center 

Operating 

Wavelength 

Width  of 
Window 

Number  of  WDM 
Channels  (4  nm  spacing) 

SI  PMMA  Plastic 

650  nm 

10  nm 

3 

SI  PMMA  Plastic 

570  nm 

35  nm 

9 

SI  PMMA  Plastic 

520  nm 

40  nm 

10 

GI  PMMA  Plastic 

650  nm 

10  nm 

3 

GI  PMMA  Plastic 

570  nm 

30  nm 

8 

GI  PMMA  Plastic 

520  nm 

0  nm 

0 

Figure  Captions 


Figure  1.  Schematic  experimental  setup  for  POF  attenuation  measurements. 

Figure  2.  Attenuation  spectra  of  (a)  SI  (1000  pm  in  diameter)  and  (b)  GI  (750  pm  in 
diameter)  PMMA  based  POF.  The  loss  at  650  nm  is  110  dB/km  for  SI  POF  and  158 
dB/km  for  GI  POF.  The  slight  increase  in  the  short  wavelength  region  for  GI  POF  loss 
spectrum  is  due  to  the  involvement  of  benzyl  benzoate. 

Figure  3.  Attenuation  spectra  of  PMMA  based  SI  (dashed  curve)  and  GI  (dotted  curve) 
POFs  and  a  theoretical  perfluorinated  POF  (solid  curve). 

Figure  4.  Schematic  of  a  5-channel  wavelength  division  multiplexing  (WDM)  system. 
Five  narrow  linewidth  lasers  transmit  the  distinct  channels  which  are  multiplexed  onto  a 
single  fiber.  After  traveling  a  long  distance  the  signals  are  demultiplexed,  regenerated, 
and  remultiplexed  back  onto  a  fiber  for  further  transmission.  At  the  terminal  end,  the 
channels  are  demultiplexed  and  individually  detected  and  processed. 

Figure  5.  Schematic  of  an  8-terminal  star  network  in  which  each  terminal,  or  user, 
transmits  a  signature  wavelength  and  can  read  all  8  wavelengths;  this  allows  each  user  to 
communicate  with  any  other  user. 

Figure  6.  Absorption  spectrum  of  Silica  glass  optical  fiber  with  transmission  windows  at 
1330  and  1550  nm. 

Figure  7.  Schematic  experimental  setup  for  time-domain  POF  bandwidth  measurement. 
The  fast  pulse  detector  is  an  optical  sampling  oscilloscope  from  Hamamatsu. 

Figure  8.  Schematic  illustration  of  a  mode  scrambler  to  achieve  equilibrium  mode 
distribution  condition  in  POFs  . 

Figure  9.  Optical  pulse  width  from  a  64  m  GI  POF  as  a  function  of  incidence  angle. 
Hollow  and  filled  circles  are  results  when  guided-modes  are  unscrambled  and  scrambled, 
respectively. 


Figure  10.  Optical  pulses  from  the  diode  (left)  and  a  12.2  m  SI  POF.  The  pulse 
broadening  is  due  to  the  intermodal  dispersion  in  the  fiber. 

Figure  11.  Optical  pulse  from  the  diode  (left)  and  a  98.8  m  GI  POF.  The  pulse 
broadening  is  due  to  the  intermodal  dispersion  in  the  fiber. 

Figure  12.  Calculated  impulse  response  function  for  an  SI  fiber.  The  spike  near  the  cut-off 
is  due  to  those  cut-off  modes  that  were  not  included  in  the  WKB  approximation. 

Figure  13.  (a).  A  five-step  stair  case  profile,  (b)  The  ratio  of  pulse  broadening  for  the  stair 
case  profile  to  that  of  the  step  index  profile. 

Figure  14.  Calculated  radial  part  of  the  electric  field  distribution  for  three  different 
guided-modes  (v,  /),  where  v  is  the  azimuthal  number  and  I  is  the  sequential  number. 

Figure  15.  Schematic  experimental  setup  of  the  refracted  near  field  measurement  of  POF 
refractive  index  profiles.  A:  chopper;  B;  polarizer;  C:  quarter-wave  plate;  D:  lOX 
microscope  objective;  E:  pinhole  (25  jim  in  diameter);  and  F;  microscope  objective.  The 
back  end  of  the  fiber  is  illuminated  with  a  He-Ne  laser  in  order  to  observe  the  front  end 
with  an  eyepiece. 

Figure  16  The  blocking  disk  assembly  of  the  refracted  near  field  measurement  technique 
of  Fig.  15. 

Figure  17.  Measured  refractive  index  profile  of  current  GI  POF  (dashed  curve)  and  its 
polynomial  fit  (solid  curve).  The  measured  curve  consists  of  three  regions,  A:  index 
matching  liquid;  B:  cladding;  and  C:  core. 

Figure  18.  Measured  refractive  index  profile  of  current  GI  POF  (dashed  curve)  and  its 
alpha  fit  (solid  curve)  with  a  =  3.67. 

Figure  19.  Calculated  impulse  response  function  based  on  the  measured  refractive  index 
profile  of  current  GI  POF  (100  m).  The  response  function  has  two  parts;  a  sharp  peak, 


corresponding  to  the  modes  confined  in  the  central-core  region,  and  a  long  tail, 
representing  the  remaining  modes  which  have  large  intermodal  dispersion. 

Figure  20.  Transfer  function  of  impulse  response  function  of  Fig.  19  based  on  the 
measured  refractive  index  profile  of  current  GI  POF(100  m).  The  bandwidth-length 
product  is  determined  to  be  450  Mhz-100  meters. 

Figure  21.  Schematic  of  the  light  paths  in  the  transverse  interferometric  method  (TIM) 
experiment. 

Figure  22.  Schematic  of  the  observed  interference  fringes  of  the  TIM  measurerrient  (a) 
before  and  (b)  after  insertion  of  the  POF  sample. 

Figure  23.  Geometry  of  the  path  traveled  by  a  ray  of  light  as  it  passes  through  the  fiber 
sample. 

Figure  24.  Diagram  of  the  sample  stage  for  the  TIM  measurement.  A  fiber  sample  is 
placed  on  a  quartz  plate  and  immersed  in  index  matching  liquid  on  the  sample  side  of  the 
stage.  The  same  index  matching  fluid  and  an  identical  quartz  plate  comprise  the 
reference  side. 

Figure  25.  Comparison  of  the  GI  POF  index  of  refraction  profile  measured  using  both  the 
transverse  interferometric  method  (solid  curve)  and  the  near  field  refracted  ray  technique 
(dashed  curve). 
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Section  VI.  Mode  Coupling  and  Equilibrium  Mode  Distribution  Conditions  in 
Plastic  Optical  Fibers 

The  most  important  property  of  a  plastic  optical  fiber  link  in  data  transmission 
applications  is  its  bandwidth,  which  defines  the  limit  of  information  carrying  capacity. 
The  bandwidth  measurement  can  be  performed  in  the  time-domain  by  measuring  its  pulse 
broadening.  *'5  However,  it  is  well  known  that  the  experimental  results  depend 
sensitively  on  the  launching  conditions,  such  as  the  incident  angle  and  numerical  aperture 
of  the  input  beam.^'^  The  reason  is  the  equilibrium  mode  distribution  (EMD)  condition 
is  not  yet  completely  achieved  in  short  fiber  waveguides  for  these  principal  optical 
characteristics  to  be  accurately  defined. 

The  EMD  condition  is  naturally  achieved  after  light  travels  a  certain  fiber  distance 
or  coupling  length  Lc,  depending  on  the  fiber  quality.  As  a  result  of  the  diffraction  of  the 

light  by  microscopic  anomalies  randomly  distributed  in  the  core,  power  transfers  from 
one  mode  to  another  until  the  power  distribution  reaches  its  equilibrium  condition.  These 
fiber  anomalies  include  extrinsic  features  such  as  micro  bends,  diameter  variations,  micro 
voids  and  cracks  in  addition  to  intrinsic  static  and  dynamic  density  fluctuations  of  the 
polymer  chains.  The  transfer  of  energy  between  modes  with  different  propagation 
velocities  tends  to  average  out  the  total  propagation  delays,  thus  reducing  the  intermodal 
dispersion  and  increasing  the  fiber  bandwidth.  For  a  fiber  long  compared  to  its  coupling 
length  Lc,  mode  coupling  is  complete  and  the  relative  mode  population  is  nearly 

independent  of  the  launching  conditions.  In  this  section,  we  summarize  the  first 
observations  of  mode  coupling  effects  in  SI  PDFs  using  a  time-domain  bandwidth 
measurement  technique,  demonstrating  that  the  coupling  length  can  be  significantly 
reduced  with  a  proper  launching  condition.  These  observations  have  also  been 
independently  verified  by  our  measurement  of  the  far  field  radiation  patterns. 

The  step  index  (SI)  fiber  used  in  our  experiment  was  ESKA  Premier  GH  4001P 
(Mitsubishi  Rayon  Co.,  Ltd.),  which  has  a  core  of  polymethylmethacrylate  (PMMA)  and 


a  cladding  of  a  fluorinated  polymer.  The  refractive  indices  of  the  core  and  cladding  are 
ni=1.492  and  n2=1.402,  respectively,  which  yields  a  theoretical  numerical  aperture  (NA) 

of  0.51.  The  fiber  diameter  is  1mm.  For  our  measurements,  the  fiber  ends  were  polished 
by  successively  using  four  dry  lapping  films  of  decreasing  size  of  5,  3,  1,  and  0.3  |im, 
respectively.  Compared  to  the  hot  plate  fiber  polishing  method,  our  technique  eliminates 
possible  end  surface  smearing  and  fiber  distortion  that  are  common  to  the  hot  plate 
method.  We  believe  our  dry  polishing  method  is  especially  suitable  for  large  core  POFs. 

Fig.  1  shows  the  experimental  setup  for  measuring  pulse  broadening  and 
transmission  bandwidth.  The  light  source  is  a  temperature  stabilized  InGaAlP  laser  diode 
(Hamamatsu  PLP-02)  capable  of  emitting  45  ps  pulses  at  660  nm  in  wavelength.  The 
repetition  rate  of  the  laser  diode  can  be  as  high  as  2MHz.  The  laser  diode  output  is  first 
collimated  by  a  GRIN  lens  and  passes  through  a  spatial  filter  to  achieve  a  decent  spatial 
beam  profile.  It  is  then  collimated  by  a  15cm  lens  and  focused  onto  the  fiber  by  one  of 
four  micro  objectives  (5X,  0.12  NA;  lOX,  0.25  NA;  20X,  0.40  NA;  40X,  0.65  NA), 
which  provide  a  source  with  controllable  launch  numerical  aperture. 

The  light  output  from  the  fiber  is  detected  by  a  Hamamatsu  sampling  optical 
oscilloscope  (OOS).  The  main  part  of  the  OOS  system  consists  of  a  sampling  streak  tube 
to  convert  the  light  incident  the  photocathode  into  photoelectrons.  These  photoelectrons 
are  sampled  and  then  directed  to  a  phosphor  screen,  where  they  are  converted  to  light. 
The  light  is  then  detected  by  a  photomultiplier  tube  (PMT).  The  data  stored  in  the 
oscilloscope  can  be  extracted  to  the  computer  for  further  data  analysis.  The  time 
resolution  of  the  OOS  system  is  better  than  10  ps.  The  pulse  broadening  (a)  is  determined 
by  the  following  equation 

cr  =  (1) 

where  x\  and  X2  are  the  full  widths  at  half  maximum  (FWHM)  of  the  light  source  and 
output  from  the  fiber,  respectively. 


Measurements  were  first  performed  to  establish  the  fiber  length  dependence  of 
the  pulse  broadening  under  different  launching  conditions.  Fig.  2  shows  the  log-log  plot 
of  the  optical  pulse  broadening  (a)  versus  the  fiber  length  (L)  for  four  different  launching 
conditions.  We  notice  that  when  the  fiber  length  is  less  than  30  m,  the  pulse  broadening 
increases  with  an  increasing  numerical  aperture,  in  agreement  with  the  understanding  that 
a  high  NA  launch  will  excite  more  guided  modes  in  the  fiber  and  thus  will  result  in  more 
pulse  broadening.  In  addition,  a  fit  to  the  following  relation 

a  (2) 

shows  that  a  =  0.99  ±  0.01  for  NA  =  0.12  and  a  =  0.97  ±  0.02  for  NA  =  0.25.  These 
results  indicate  that  for  under-filled  launching  cases,  the  pulse  broadening  is  nearly  a 
linear  function  of  the  fiber  length.  A  slight  decrease  in  the  exponent  a  with  an  increased 
NA  suggests  that  the  mode  coupling  has  already  started,  even  in  an  SI  POF  as  short  as 
20m.  More  importantly,  when  the  light  is  coupled  into  the  fiber  using  a  high  NA 
microscope  objective,  the  pulse  broadening  clearly  shows  two  distinct  regions.  In  the 
first  region,  a  fit  to  Eq.  (2)  yield  a  =  0.97  ±  0.03  for  NA  =  0.40  and  a  =  0.96  ±  0.10  for 
NA  =  0.65,  indicating  the  mode  coupling  is  far  from  being  complete.  In  the  second 
region,  the  pulse  broadening  slows  down  significantly,  with  the  fit  producing  a  =  0.54  ± 
0.01  for  NA  =  0.40  and  a  =  0.57  ±0.01  for  NA  =  0.65. 

A  diffusion  theory  of  mode  coupling  in  a  fiber  has  been  established  by  Gloge^. 
The  theory  predicted  that  when  the  mode  coupling  is  complete  and  equilibrium  mode 
distribution  (EMD)  is  established,  the  optical  pulse  broadening  changes  from  being 
linearly  proportional  to  the  fiber  length  (a  =  1)  to  varying  as  the  square  root  of  the  length 
(a  =  0.5).  Our  results  clearly  shows  that  the  agreement  between  the  diffusion  theory  and 
experiment  is  quite  satisfactory.  In  addition,  our  measurement  indicates  that  the  coupling 
length  is  around  15  to  20  m.  The  much  shorter  coupling  length  compared  to  glass  optical 
fibers  which  have  a  typical  coupling  length  of  several  kilometers^  shows  that  the  mode 
coupling  effects  are  much  more  significant  inside  POFs.  This  coupling  length  can  further 


be  shortened  if  the  light  injected  into  the  fiber  has  a  spatial  profile  close  to  the  one  inside 
the  fiber  when  the  EMD  condition  has  been  achieved. 

We  also  carried  out  far  field  pattern  measurements  to  study  the  mode  coupling 
effect  under  different  launching  conditions.  Measurements  were  performed  on  a  30  m- 
long  SI  POP  using  a  CW  LED  at  650nm.  After  being  modulated  by  a  chopper,  the  light  is 
focused  onto  the  input  end  of  the  fiber  mounted  on  a  rotation  stage  which  allows  the 
launch  angle  to  be  varied  with  respect  to  the  fiber  axis.  The  output  end  of  the  fiber  was 
mounted  at  the  center  of  another  rotation  stage.  Angular  scans  of  the  far  field  intensity 
were  made  using  a  silicon  photodiode  with  a  small  active  area  of  Imm^,  which  was  set 
on  an  arm  that  can  rotate  around  the  output  fiber  end.  The  photodiode  and  the  fiber  end 
were  15  cm  apart,  so  the  angular  resolution  was  better  than  half  a  degree.  The  electrical 
signal  from  the  photodiode  was  sent  to  a  lock-in  amplifier  referenced  by  the  chopper. 

Fig.  3  and  4  show  angular  scans  of  the  far  field  intensity  under  a  variety  of  launch 
conditions.  For  each  experiment,  light  was  launched  into  the  fiber  at  angle  with  respect  to 
the  fiber  axis  of  from  0  to  25  degrees  in  5  degree  increments.  In  Fig.  3,  the  light  was 
launched  into  the  fiber  using  a  low  numerical  aperture  (NA  =  0.12)  microscope 
objective.  We  can  clearly  see  that  the  far  field  pattern  changes  from  being  a  disk  to  a  ring 
when  the  incidence  angle  exceeds  20°,  suggesting  that  the  only  a  group  of  high  order 
modes  is  excited,  in  agreement  with  our  previous  finding  in  the  time  domain  bandwidth 
measurement  that  the  mode  coupling  is  far  from  being  complete.  In  contrast,  as  can  be 
seen  in  Fig.  4,  when  the  light  was  launched  into  the  fiber  with  a  high  NA  (0.40) 
microscope  objective,  the  output  distribution  is  essentially  independent  of  the  launch 
condition  and  consists  of  a  disk,  regardless  of  the  incident  angle  provided  it  is  within  the 
numerical  aperture  of  the  fiber,  suggesting  that  the  mode  coupling  is  complete  and  the 
EMD  condition  is  achieved.  These  results  firmly  confirmed  our  bandwidth  observations. 


According  to  the  WKB  methodic,  ^hich  is  a  good  approximation  as  the  number  of 

guided  modes  in  the  SI  POF  is  several  millions,  the  modified  mode  number  m  is  related 

to  the  incidence  angle  9  by  the  following  equation^ 2 
m  _  sin  6 

M  sin  (3) 

where  M  is  the  total  number  of  guided  modes  in  the  fiber  and  0c  is  the  acceptance  angle 
of  the  fiber  determined  by  its  numerical  aperture.  Therefore,  the  transmission  of  a  fiber 
as  a  function  of  the  incidence  angle  is  related  to  its  differential  attenuation  coefficient 
provided  the  mode  coupling  is  weak.  In  Fig.  3,  the  area  under  a  curve  is  proportional  to 
the  total  light  intensity  coming  from  the  output  end.  We  can  therefore  conclude  that  high 
order  modes  exhibit  higher  attenuation  in  the  fiber,  as  indicated  from  Fig.  3.  The  slight 
decrease  in  light  intensity  with  an  increasing  incidence  angle  in  Fig.  4  is  due  to  the  fact 
that  modes  in  the  fiber  requires  about  15  to  20  m  before  they  are  fully  coupled. 

In  summary,  we  have  carried  out  pulse  broadening  measurements  of  SI  POFs  as  a 
function  of  fiber  length  under  different  launching  conditions,  with  the  results  showing 
that  the  HMD  condition  can  readily  be  achieved  in  an  SI  POF  much  shorter  than  in  glass 
optical  fibers.  The  coupling  length  for  the  POF  (NA  =  0.51)  was  determined  to  be  around 
15  to  20m  for  a  launching  NA  of  0.40.  There  is  a  power  law  dependence  between  the 
pulse  broadening  and  the  fiber  length,  with  the  measured  exponent  close  to  1 .0  before  the 
HMD  condition  is  achieved  and  0.5  after  the  HMD  condition  is  established.  The 
observation  of  the  HMD  condition  was  also  independently  verified  using  the  far  field 
pattern  technique. 
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Figure  Captions 

Fig.  1.  Schematic  experimental  setup  for  time-domain  POF  bandwidth  measurement. 
M.O.rmicro  objective;  OOS:  Hamamatsu  optical  sampling  oscilloscope. 

Fig.  2.  Log-log  plot  of  output  pulse  width  as  a  function  of  fiber  length  with  launch 
numerical  aperture  as  the  parameter. 

Fig.  3.  Far-field  radiation  patterns  from  a  30m-long  POF  under  low  NA  launch 
conditions,  launch  NA=0.12. 

Fig.  4.  Far-field  radiation  patterns  from  a  30m-long  POF  under  high  NA  .launch 
conditions,  launch  NA=0.4. 
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Section  VII.  Origin  of  High  Bandwidth  Performance  of  Graded-Index  Plastic 
Optical  Fibers 

The  main  limitation  on  the  bandwidth  of  multimode  fibers’’^  is  due  to  modal 
dispersion  in  which  different  optical  modes  propagate  at  different  group  velocities.  There 
are  essentially  two  principal  means  to  compensate  for  the  dispersion,  namely  by  forming 
a  parabola-like  graded  refractive  index  profile  across  the  fiber  core  and  by  optical  mode 
coupling  in  which  energy  is  exchanged  between  the  modes  by  way  of  random 
perturbations  along  the  propagation  direction  of  the  fiber.  Up  to  the  present  time,  it  has 
been  thought  that  the  increased  bandwidth  performance  of  GI POF  arises  strictly  from  the 
G1  profile.  In  this  work,  we  show  that  the  index  profile  of  the  current  GI  POF  is  not 
parabolic  and  that  the  bandwidth  performance  is  largely  determined  by  strong  mode 
coupling. 

The  GI  POF  samples  were  fabricated  by  standard  melt-draw  near  180  -  200  °C  at 
a  pulling  rate  of  1  m  per  3  min  from  molecularly  doped  poly(methyl  methacrylate) 
(PMMA)  GI  preforms  made  as  previously  described^  A  PMMA  tube  was  formed  by 
heating  a  rotating  glass  tube  at  70  -  80  °C  containing  highly  purified  methyl  methacrylate 
(MMA)  monomer,  benzyl  peroxide  initiator,  and  n-butyl  mercaptan  chain  transfer  agent 
and  heat  treated  near  100  °C  at  reduced  pressure  overnight.  The  polymerized  PMMA 
tube  was  filled  with  mixture  of  MMA  monomer  and  approximately  20  %  wt  benzyl 
benzoate  molecular  guest  along  with  suitable  small  amounts  of  initiator  and  chain  transfer 
agent  and  then  heat-treated  at  about  70  -  90  °C.  All  the  fabricated  fibers  possessed 
optical  losses  less  than  175  dB/km  at  650  nm.  The  core  of  each  fiber  is  approximately 
500  ^im  in  diameter  and  consists  of  a  mixture  of  PMMA  and  benzyl  benzoate  which  has  a 
higher  refractive  index  than  that  of  PMMA.  The  cladding  layer  is  125  pm  thick  and 
comprised  of  pure  PMMA  to  yield  a  numerical  aperture  (NA)  of  0.25  for  each  fiber. 

The  bandwidth  of  the  each  GI  POF  samples  was  measured  in  the  time  domain 
using  a  temperature  stabilized  InGaAlP  laser  diode  that  produces  45  ps  pulses  at  660  nm 
with  a  repetition  rate  of  2  MHz.  The  laser  diode  output  is  collimated  by  a  GRIN  lens, 
spatially  filtered,  and  then  collimated  by  a  15  cm  focal  length  lens  to  be  focused  onto  the 


POF  end  by  a  microscope  objective  lens  (20X,  0.40  NA).  The  fiber  output  is  detected  by 
a  sampling  optical  oscilloscope.  Figure  1  shows  the  input  and  output  optical  pulses  for  a 
GI  POF  98.8  m  in  length.  From  fast  Fourier  transforms  of  the  pulses,  the  measured 
bandwidth  is  determined  to  be  3.0  GHz  for  100m  length  fiber. 

For  the  refractive  index  profile  measurements,  we  implemented  two  independent 
optical  techniques;  the  refracted  near  field'*’^  and  transverse  interferometric  methods^.  The 
results  of  both  methods  were  consistently  in  agreement.  In  the  latter  method  which 
utilizes  an  interference  microscope,  the  displacements  of  interference  fringes  generated  by 
a  fiber  sample  are  recorded  by  a  CCD  camera  and  used  to  calculate  the  refractive  index 
profile  of  the  fiber.  The  measured  profile  for  a  typical  GI  POF  sample  is  plotted  in  Figure 
2  as  the  dashed  curve.  The  distinguishing  feature  of  the  measured  profile  is  that  it  is  not 
parabolic  as  one  would  expect  for  a  true  GI  fiber.  The  central  part  of  the  profile,  close  to 
the  core  center,  is  closely  parabolic,  but  the  outer  region  connects  steeply  to  the  cladding, 
reminiscent  of  low  bandwidth  multimode  step  index  (SI)  fiber.  Additional  GI  POF 
samples  fabricated  independently  following  the  same  methods  all  showed  similar  index 
profiles.  This  demonstrates  that  the  characteristic  features  of  the  refractive  index  profile 
are  intrinsic  to  the  preform  fabrication  process. 

Initially,  we  attempted  to  fit  the  measured  data  with  an  alpha  profile  and  even  a 
single  sixth-order  polynomial,  and  found  that  these  functions  are  not  able  to  adequately 
represent  important  features  throughout  the  entire  core  region.  We  find  that  two 
polynomials,  corresponding  to  the  central-core  (r  <  aO  and  core-cladding  {a\  <  r  <  a) 
regions,  are  necessary  to  properly  fit  the  measured  profile.  Accordingly,  we  represent  the 
profile  with  the  function 


«(r)  = 


n. 


n. 


-Sn^Q(r^  a)^  —dn^^{r  I  a)‘^ 
-Sn2Q(rl a)^  -  Sn2^(r  I a)^ 


for  r  <  a, , 
for  a^<r<a, 
for  r>  a. 


(31) 


With  the  requirement  that  the  polynomials  connect  smoothly  at  r  =  a\,  the  final 
parameters  are  =  1.51 19,  bnw  =  0.0077,  5nii  =  0.0088,  =  1.5086, 8n2o  =  -  0.0073, 


5/121  =  0.0256,  /ic  =  1.4910,  fli  =  161.4  |im,  and  a  =  257  p.m.  The  polynomial  fit  is  plotted 
in  Figure  2  as  the  solid  curve.  Figure  3  shows  the  measured  profile  along  with  its  best 
alpha  profile  fit  having  an  exponent  of  a  =  3.64.  The  alpha  profile  does  not  agree  with  the 
data  in  the  central  core  region. 

Determination  of  the  bandwidth  from  the  measured  refractive  index  profile 
requires  an  estimate  of  the  modal  power  distribution  inside  the  fiber.  Following  the 
model  typically  used  for  graded-index  glass  fibers’,  the  modal  power  distribution  p(m)  for 
our  calculations  is  taken  to  be 


p(m)  = 


fl -(/n/M,) 

[0 


for  m  < 
for  m  > 


(32) 


where  m  is  the  number  of  modes  above  the  propagation  constant,  (3.  Me  represents  the 
maximum  number  of  modes  that  can  propagate  in  the  fiber.  Here,  Me  is  set  to  0.95M, 
where  M  is  the  total  number  of  guided  modes.  The  pulse  broadening  (impulse  response 
function,  h{x))  is  then  calculated,  taking  into  account  intermodal  dispersion,  using  the 
following  equations  from  the  WKB  approximation*: 


m(fi,k)  =  j'^[n\r)k^-j3^]rdr. 


(33) 
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where  m  is  again  the  number  of  modes  above  the  propagation  constant,  p,  k  is  the  wave 
number,  r^  is  determined  by  ni^r^ )  =  pi  k  ,  and  x  is  the  group  delay  time.  The  above 
equations  assume  mode  coupling  is  not  present. 

In  principle,  an  ideal  parabolic  index  profile  fully  compensates  for  the  longer 
distances  traveled  by  rays  of  the  outer  modes  so  that  all  guided  modes  propagate  at  the 


same  speed,  essentially  eliminating  intermodal  dispersion.  Since  the  refractive  index 
profile  of  the  present  GI  POP  most  closely  resembles  an  alpha  profile  greater  than  the 
ideal  parabolic,  it  under-compensates  for  the  outer  modes  near  the  cladding.  In  other 
words,  the  local  refractive  index  approaching  the  cladding  is  too  high  such  that  rays  in  the 
outer  modes  travel  at  slower  speeds  than  those  toward  the  fiber  center,  resulting  in 
increased  intermodal  dispersion.  As  such,  material  dispersion  is  not  dominant  since  it  is 
small  in  comparison  to  the  intermodal  dispersion.  In  particular,  we  find  for  alpha  profiles 
with  a  =  3.64  and  source  spectral  widths  of  less  than  5  nm,  as  realized  in  our 
experiments,  material  dispersion  increases  the  bandwidth  by  only  12.5%. 

The  Fourier  transform  of  the  impulse  response  function  h(T)  yields  the  bandwidth, 
which  for  the  current  GI  POP  is  determined  to  be  0.38  GHz  for  100  m,  with  an  additional 
0.05  GHz  from  material  dispersion.  Other  modal  power  distribution  models  have  been 
used,  including  distributions  with  different  values  of  Me  as  well  as  an  equal  modal  power 
distribution.  They  all  result  in  calculated  bandwidths  much  smaller  than  the  measured 
bandwidth  of  3.0  GHz  for  100m.  We  ascribe  the  bandwidth  difference  between  the 
measured  value  of  3.0  GHz  and  the  profile-determined  value  of  0.43  GHz  as  being  due  to 
strong  mode  coupling  that  effectively  increases  the  bandwidth. 

A  direct  test  for  the  presence  of  strong  mode  coupling  is  to  measure  the 
bandwidth,  or  equivalently  the  inverse  pulse  broadening  dependence  (a)  on  fiber  length 
(L)  given  as  a  L^.  In  a  multimode  fiber,  mode  coupling  is  the  switching  back  and  forth 
of  energy  packets  between  different  modes.  If  there  is  no  mode  coupling,  the  difference 
in  arrival  times  between  the  energy  packets  in  the  fast  mode  and  those  in  the  slow  mode  is 
proportional  to  the  fiber  length.  However,  if  the  energy  packets  randomly  switch  back 
and  forth  between  the  fast  and  slow  modes,  they  arrive  at  the  end  of  the  fiber  closer 
together  in  time.  In  fact,  it  can  be  shown  that  the  width  of  the  distribution  of  energy 
packet  arrival  times  is  proportional  to  the  square  root  of  the  fiber  length.  A  detailed  mode 
coupling  theory  shows  that,  once  an  equilibrium  mode  distribution  is  established,  pulse 
broadening  increases  as  the  square  root  (P  =  0.5)  of  length  rather  than  as  the  usual  linear 
(P  =  1)  dependence.  The  length  at  which  the  pulse  broadening  dependence  changes  from 


P  =  1  to  P  =  0.5  is  called  the  coupling  length,  L,..  We  performed  these  pulse  broadening 
measurements  on  GI  POFs,  using  the  same  setup  as  in  the  time  domain  bandwidth 

measurements.  The  pulse  broadening  is  defined  by  the  relation  g  =  ^[T2  ~  '^\  >  where  X| 

and  X2  are  the  full  widths  at  half  maximum  (FWHM)  of  the  input  and  output  pulses, 
respectively.  A  log-log  plot  of  the  experimentally  determined  pulse  broadening  as  a 
function  of  fiber  length  is  given  in  Figure  4  along  with  a  least  squares  fit  (solid  line).  The 
measured  data  indeed  exhibit  strong  mode  coupling  behavior  with  the  experimental  value 
p  =  0.6,  close  to  the  ideal  strong  coupling  limit  of  p  =  0.5  and  far  from  the  weak  coupling 
linear  dependence  of  P  =  1. 

Our  results  also  demonstrate  that  the  coupling  length  of  current  GI  POFs  is  shorter 
than  the  experimental  lower  limit  of  30  m.  Furthermore,  it  can  be  shown  that 


where  and  Bko  are  the  bandwidths  in  the  presence  and  absence  of  strong  mode 
coupling,  respectively.  Using  our  measured  and  profile  determined  values  of  3.0  and  0.43 
GHz,  respectively,  for  100  m,  the  value  of  Lc  is  determined  to  be  2.1  m,  which  is 
consistent  with  experiment.  In  contrast  to  kilometer  long  coupling  lengths  characteristic 
of  GOF,  such  short  observed  Lc  reflects  intrinsic  density  and  concentration  fluctuations  in 
addition  to  the  usual  extrinsic  perturbations  such  as  microbends.  The  presence  of  the 
material  fluctuations  is  consistent  with  the  low  glass  transition  temperatures  observed  in 
GI  POFs. 

In  conclusion,  we  have  made  the  first  observations  in  graded  index  plastic  optical 
fibers  (GI  POFs)  of  strong  mode  coupling  which  largely  determines  the  measured  high 
bandwidth  performance.  We  observed  the  characteristic  square  root-like  fiber  length 
dependence  of  the  optical  pulse  broadening.  We  have  further  found  that  the  measured 
refractive  index  profile  is  not  parabolic  and  is  best  described  by  a  dual  polynomial  in 
which  the  central-core  and  outer  core-cladding  regions  are  markedly  different.  Near  the 


central  core  the  profile  is  closely  parabolic,  but  the  outer  region  drops  steeply  to  the 
cladding,  similar  to  a  step  index  fiber.  Analysis  of  the  index  profile  reveals  that  strong 
mode  coupling  increases  the  GI  POP  bandwidth  from  its  profile-determined  value  of  0.43 
GHz  to  its  measured  value  of  3.0  GHz  for  100  m.  Large  density  and  concentration 
fluctuations  natural  to  plastic  fiber,  in  addition  to  extrinsic  perturbations,  are  responsible 
for  the  relatively  short  coupling  lengths  estimated  to  be  around  2m.  Lastly,  recent 
independent  studies  of  differential  mode  delay  and  attenuation  completed  after  the 
DARPA  grant  period  largely  confirm  the  refractive  index  profile  and  bandwidth  results 
for  GI  POP.  The  results  of  these  new  studies  will  be  published  during  1998. 
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Figure  Captions 

Figure  1.  Measured  input  optical  pulse  (dashed  curve)  and  output  pulse  (solid  curve) 
from  a  98.8  m  GI  POF.  The  intensity  of  both  curves  has  been  normalized.  The  full 
widths  at  half  maximum  are  45  and  145  ps  for  the  input  and  output  pulses,  respectively.  A 
fast  Fourier  transform  is  used  to  calculate  the  bandwidth  of  3.Q±0.3  GHz  for  100  m. 

Figure  2.  Refractive  index  profile  for  a  GI  POF  measured  using  a  transverse 
interferometric  method  (dashed  curve)  and  its  dual  polynomial  fit  (solid  curve).  The 
wavelength  at  which  the  measurement  is  performed  is  546  nm.  An  independent 
measurement  using  a  refracted  near-field  technique  at  650  nm  yielded  nearly  identical 
results. 

Figure  3.  Measured  refractive  index  profile  for  a  GI  POF  (dashed  curve)  and  its  alpha  fit 

(solid  curve)  with  a  =  3.64.  An  alpha  profile  is  defined  as  n(r)  =  ~('^/  ’ 

where  n/  and  tic  are  the  refractive  indices  at  the  fiber  core  center  and  at  the  cladding, 
respectively,  and  a  is  the  fiber  core  radius. 

Figure  4.  Log-log  plot  of  pulse  broadening  as  a  function  of  fiber  length  for  a  typical  GI 
POF  (points)  and  its  fit  to  <J<=<  ,  with  ^  =  0.6  ±  0.1 . 
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